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Magnetic behaviour of a mixed spin-1/2 and spin-1 Ising model on the diced 
lattice is studied by the use of an exact star-triangle mapping transformation. It 
is found that the uniaxial as well as biaxial single-ion anisotropy acting on the 
spin-1 sites may potentially cause a reentrant transition with two consecutive 
critical points. Contrary to this, the effect of next-nearest-neighbour interaction 
between the spin-1/2 sites possibly leads to a reentrant transition with three 
critical temperatures in addition to the one with two critical points only. The 
shape of the total magnetization versus temperature dependence is particularly 
investigated for the case of ferrimagnetically ordered system. 
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1. Introduction 

Over the last six decades, much effort has been devoted to determine a criticahty 
and other statistical properties of various lattice-statistical models, which would 
enable a deeper understanding of order-disorder phenomena in solids. The planar 
Ising model takes a prominent place in the equilibrium statistical physics from this 
point of view as it represents a rare example of exactly solvable lattice-statistical 
model since Onsager's pioneering work [Q]. It is noteworthy that an exact tractability 
of two-dimensional Ising systems [ |2] is being of immense importance because of 
providing a convincing evidence for many controversial results predicted in the phase 
transition theory and moreover, the exactly soluble models provide a good testing 
ground for approximative theories. It should be emphasized, however, that a precise 
treatment of planar Ising models is usually connected with the usage of sophisticated 
mathematical methods, which consequently lead to considerable difficulties when 
applying them to more complicated but simultaneously more realistic systems. 

*This work was financially supported under the grants VEGA 1/2009/05 and APVT 20-005204. 
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The Ising systems consisting of mixed spins of different magnitudes, so-called 
mixed-spin Ising models, are being among the most interesting extensions of the 
standard spin-1/2 Ising model which are currently of great research interest. These 
models enjoy a considerable attention mainly due to much richer critical behaviour 
they display compared with their single-spin counterparts. Magnetic properties of 
the mixed-spin Ising models can essentially be modified, for instance, by a presence 
of single-ion anisotropy acting on the S > 1 spins. Indeed, this effect induced by a 
crystal field of ligands may potentially cause a tricritical phenomenon, a change of 
the character of phase transition from a second-order to a first-order one. Another 
aspect which started to attract an appreciable interest towards the mixed-spin Ising 
models can be related to a theoretical modeling of magnetic structures suitable for 
describing a ferrimagnetism of certain class of insulating materials. Accordingly, the 
ferrimagnetic mixed-spin Ising models are very interesting also from the experimen- 
tal viewpoint, first of all in connection with many possible technological applications 
of ferrimagnets to practice (e.g. thermomagnetic recording). 

Despite the intensive studies, there are only few examples of exactly solvable 
mixed-spin models, yet. Using the generalized form of decoration-iteration and star- 
triangle transformations, the mixed spin-1/2 and spin-S" {S > 1) Ising models on the 
honeycomb, diced and decorated honeycomb lattices [0], as well as on the decorated 
square, triangular and diced lattices [E] were exactly been treated long ago. Notice 
that these mapping transformations have been later on further generalized in order to 
account also for the single-ion anisotropy effect. The influence of uniaxial single-ion 
anisotropy on magnetic properties of the mixed-spin systems have been precisely 
investigated on the honeycomb lattice [ E] and on some decorated planar lattices 
[ E] . In addition, more general form of the star-triangle transformation regarding 
the biaxial single-ion anisotropy have only recently been employed to obtain exact 
results for the mixed-spin Ising model on the honeycomb lattice [|7]. To the best of 
our knowledge, these are the only mixed-spin systems with generally known exact 
solutions except the ones on the Bethe (Cayley tree) lattices, which can be studied 
within a discrete non-linear map [|H] and/or an approach based on exact recursion 
equations [IH]. Among the remarkable models for which a precise solution is restricted 
to a certain subspace of interaction parameters only, one should further mention the 
mixed-spin Ising models on the union jack (centered square) lattice that can be 
mapped onto the eight- vertex model [llUj. 

The purpose of this work is to provide an exact formulation of the mixed spin- 
1/2 and spin-1 Ising model on the diced lattice and to establish accurate phase 
diagrams of this system in dependence on a strength of the uniaxial and biaxial 
single-ion anisotropies as well as on a strength of the next-nearest-neighbour in- 
teraction. Exact results of this system are obtained by applying the generalized 
star-triangle transformation constituting an exact correspondence with the effective 
spin-1/2 Ising model on the triangular lattice. Last, the effect of single- ion anisotropy 
terms and the next-nearest-neighbour interaction on the temperature dependence 
of the resultant magnetization will be clarified for the case of ferrimagnetic system. 

The organization of this paper is as follows. In Section 2, a detailed formula- 
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Oo-=±l/2 •5'=±1,0 



Figure 1. Diagrammatic representation of the diced lattice consisting of the mixed 
spin-1/2 (empty circles) and spin-1 (solid circles) sites, respectively. The solid 
(broken) lines depict the nearest-neighbour (next-nearest-neighbour) interactions. 
By the use of star-triangle transformation, the mixed-spin diced lattice is mapped 
onto a simple spin-1/2 triangular lattice. 



tion of the model system is presented and subsequently, exact expressions for the 
magnetization and critical temperatures are derived by establishing an exact corre- 
spondence with the effective spin-1/2 Ising model on the triangular lattice. The most 
interesting numerical results are presented and particularly discussed in Section 3. 
Finally, some concluding remarks are given in Section 4. 



2. Model and its exact solution 

Let us begin by considering the mixed spin-1/2 and spin-1 Ising model on the 
diced lattice schematically illustrated on the left-hand side of figure 1. The mixed- 
spin diced lattice consists of two interpenetrating sub-lattices A and B which are 
formed by the spin-1/2 (empty circles) and spin-1 (solid circles) atoms, respectively. 
The Ising Hamiltonian defined upon the underlying lattice reads: 

6iV 3N 2N 2N 

n,.ce, = -J E s!a]-j' E aiaf - DY^isn' - EY.[is:r - {sm, (i) 

(i,j)cC {k,l)cK. i=l i=l 

where (j| and 5," (a = x, y, z) stand for the spatial components of spin-1/2 and spin- 
1 operators, respectively, J denotes the exchange interaction between the nearest- 
neighbouring spin pairs and J' labels the interaction between the next-nearest- 
neighbouring spin-1/2 atoms. Finally, the parameters D and E measure a strength 
of the uniaxial and biaxial single- ion anisotropics acting on the spin-1 sites and N 
denotes the total number of the spin-1/2 sites. 

It can be clearly seen from figure 1 that the considered diced lattice belongs to 
loose-packed lattices, i.e. it can be decomposed into two interpenetrating sub-lattices 
A and B in such a way that all nearest neighbours of the spin-1 sites are being the 
spin-1/2 sites and vice versa. Consequently, a trace over spin degrees of freedom of 
all spin-1 atoms can be performed before summing over all available configurations 
of the spin-1/2 atoms. It is therefore very convenient to write the total Hamiltonian 
(H} as a sum of the commuting site Hamiltonians Hdiced = where each site 
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Hamiltonian Tii involves all interaction terms associated with one appropriate spin-1 
site residing ith lattice point of the sub-lattice B: 

n,= - StJial, + at, + a^3) - {SffD - [{S^f - {S^E 

- J'{al^l, + ayrl, + al,al^l2 (2) 

The factor | emerging in the last expression arises in order to avoid a double count- 
ing of each next-nearest-neighbour interaction when summing over all spin-1 sites. 
On account of the commutability between the site Hamiltonians (j21), the partition 
function of the mixed-spin diced lattice can be partially factorized to: 

2N 

Z diced = Tr{^} n Tr exp ( -(iUi ) , (3) 

4 = 1 

where /3 = 1/(^3^), is being Boltzmann's constant, T stands for the absolute 
temperature, the symbol Tr5. means a trace over the spin degrees of freedom of ith 
spin-1 atom of the sub-lattice B and the symbol Trjo-} stands for a trace over all 
possible spin configurations of the spin- 1/2 atoms on the sub-lattice A. 

After performing elementary calculations (see for instance reference [E]); one 
readily finds an analytical expression for Tr^- exp(— Moreover, an explicit form 
of the relevant trace immediately implies a possibility of performing the familiar 
star-triangle mapping transformation [|3j: 

IY5,exp(-/3H,) = exp[/3J'(afia^2 + ^r2^f3 + ^r3^ri)/2)] 

{l + 2 exp(/?D) cosh[/3v' J2(afi + af2 + afs)^ + E''\ } 

= C exp[/?i?(a^,a^2 + ^l2^lz + ^-a^fJ/S] (4) 

which establishes an exact correspondence between the mixed-spin Ising model on 
the diced lattice and its equivalent spin-1/2 Ising model on the triangular lattice 
provided that the transformation (ji} is applied to all spin-1 sites (see figure 1). Notice 
that the unknown mapping parameters C and R are unambiguously given directly 
by the mapping transformation (j3)), which must hold for all available configurations 
of crfi, (jf2, o"^^3 spins. According to this, one easily finds: 

= ViVi, I3R = pJ' + 2 ln{Vi/V2), (5) 

with the expressions Vi and V2 defined as: 

Vi = 1 + 2 exp{(3D) cosh(/3y'(3 J/2)2 + E^) , (6) 
V2 = 1 + 2 expipD) cosh(/?y^( J/2)2 + E^) . (7) 

Now, let us substitute the mapping relation (jD) into the formula gained for 
the partition function of the mixed-spin diced lattice. After straightforward rear- 
rangement, one easily obtains an exact relationship between the partition function 
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Zdiced of the mixed-spin Ising model on the diced lattice and respectively, the par- 
tition function Ztnang of its corresponding spin-1/2 Ising model on the triangular 
lattice: 

Zdiced{l3-i J, J' -I D, E) = C"^^ ZtriangiP, R)- (§) 

With regard to this, the mixed-spin diced lattice is mapped onto the spin-1/2 Ising 
model on the triangular lattice exactly solved several years ago [ITT]-[IT^. 

The mapping relation (jH)) between the partition functions Zdiced and Ztriang rep- 
resents a central result of our calculation, in fact, this relation can be utilized for 
establishing similar mapping relations also for other important quantities such as 
Gibbs free energy, internal energy, magnetization, correlation functions, specific heat, 
etc. When combining the formula (jHl) with commonly used mapping theorems [I14j. 
one readily proves a validity of following relations: 

= {^i ) diced {^i ) triang — ^triang y (9) 
tA = {O'i O'j al) diced = {O'i O'j a^) triang = ttriang, (10) 

rriB = {S^)d^ced = SmA[F{3) + F{l)]/2 + 2tA[Fi3)~3F{l)], (11) 

where the symbols {■■■) diced and (...) triang denote standard canonical average per- 
formed over the Ising models on the mixed spin-1/2 and spin-1 diced lattice and on 
the spin-1/2 triangular lattice, respectively, and the function F{x) stands for: 



p ^ Jx/2 2smhi/3^iJx/2y + E^) 

'{Jx/2y + E^ exp{-(3D) + 2 cosh(/3./( Ja;/2)2 + E^) ' 



It should be pointed out that the set of equations (Pl- (fTT|l enables straightforward 
calculation of both sub-lattice magnetization niA and tub of the mixed-spin system 
on the diced lattice. According to equation the sub-lattice magnetization 
directly equals to the corresponding magnetization rritriang of the triangular lattice, 
which is being unambiguously given by the mapping relation (0). Note that the 
exact result for magnetization of the spin-1/2 Ising model on the triangular lattice 
was firstly derived by Potts in his short note [[T2|j: 



1 V W 



rrif 



riang 



1 - 



;i + 3a;)(l-x) 



(13) 



where x = exp(— /5i?). On the other hand, an exact knowledge of the triplet order 
parameter ttriang is being required in order to compute the sub-lattice magnetization 
in addition to the single-site magnetization rritriang. Exact results of Baxter and 
Choy [T^ for the anisotropic triangular Ising model include, as specialization, the 
following exact solution for the triplet correlation of the isotropic triangular lattice: 
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In view of the results (fT!^ and (fT^ . an exact solution of both sub-lattice magneti- 
zation of the mixed-spin system is thus formally completed. 

Finally, let us derive an analytical condition that determines a criticality of the 
mixed-spin diced lattice. It can easily be understood that both sub-lattice magnetiza- 
tion rriA and vanish if the magnetization rritriang of the corresponding triangular 
lattice tends to zero. It is therefore sufficient to substitute an exact critical point of 
the isotropic triangular lattice [|^ into the mapping relation (0) in order to obtain 
an accurate condition for the criticality of the mixed-spin Ising model: 

a J + 2 In "^^^ V / = In 3, (15) 

I l + 2exp(/3,D)cosh(/3,y^(J/2)2 + E2j J 

where /3c = \l{k^T^ and T^. denotes the critical temperature. 



3. Numerical results and discussion 



Before proceeding to a discussion of the most interesting numerical results it is 
worthwhile to mention that the results obtained in the preceding section hold for 
the ferromagnetic (J > 0) as well as antiferromagnetic (J < 0) nearest-neighbour 
interaction J. For simplicity, we shall restrict ourselves in what follows to the lat- 
ter case corresponding to the ferrimagnetically ordered system. It is noteworthy, 
however, that the critical frontiers of the ferromagnetic model are congruent with 
those displayed below for the ferrimagnetic model as a result of an invariance of 
the critical condition (fTH|) with respect to the transformation J — >■ —J. It should 
be also mentioned that the critical condition ()15p does not depend on a sign of the 
parameter hence, one also may consider without loss of generality E' > 0. 



3.1. Ground-state phase diagrams 

Let us firstly focus on the ground-state phase diagrams of the system under 
investigation. As the easy-plane uniaxial single-ion anisotropy (D < 0) reinforces, 
a first-order phase transition from the ferrimagnetically ordered phase towards the 
disordered phase should be expected to occur owing to energetic favoring of a non- 
magnetic state S" = of the spin-1 atoms. It is quite interesting to ascertain, however, 
that the biaxial single-ion anisotropy stabilizes the ferrimagnetically ordered ground 
state in that it enhances a critical value of the uniaxial single-ion anisotropy at which 
both ferrimagnetic as well as disordered phases coexist together: 



1^1 2^1 -{E/Dy 

This remarkable feature is in an obvious contradiction with naive expectations, 
since the biaxial term E entails uprise of local quantum fiuctuations that manifest 
themselves in a quantum reduction of the sub-lattice magnetization m^. Actually, 
arbitrary but non-zero biaxial single-ion anisotropy destroys a perfect ferrimagnetic 
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order even at zero temperature and the occurring quantum ferrimagnetic phase can 
be characterized through its sub-lattice order parameters given by: 

yriA = --, mB = , = (17) 

2 ^9 + A{E/\J\y 

As far as the influence of next-nearest-neighbour interaction J' is concerned, one 
should vigorously distinguish between the effect of ferromagnetic (J' > 0) and anti- 
ferromagnetic (J' < 0) interaction. The former interaction aligns ferromagnetically 
the next-nearest-neighbouring spin- 1/2 atoms and thus, the spin system remains 
long-range-ordered independently of a strength of the uniaxial single-ion anisotropy 
that acts purely on the spin-1 sites. On the other hand, the latter interaction causes 
a spin frustration between the next-nearest-neighbours and in the consequence of 
that, it acts in accompliance with the uniaxial single-ion anisotropy in view of de- 
stroying the ferrimagnetic long-range order. This synergetic effect can also be viewed 
from a modified version of the critical condition (jl6|) when switching on the next- 
nearest-neighbour interaction J': 

^ = , - 4^ (for J' < only). (18) 

1^1 2Jl-{E/Dy 2|J| ^ ^ ^ 



3.2. Finite-temperature phase diagrams 



Now, let us turn to a detailed discussion dealing with the finite temperature phase 
diagrams. For this purpose, figure 2 shows the variation of reduced critical tempera- 
ture with the uniaxial single-ion anisotropy for several fixed values of the ratio E/D 
and J' /\J\ = 0.0. The quantum ferrimagnetic phase can be located inside of the de- 
picted phase boundaries, while above them the usual paramagnetic phase becomes 

stable. There are two inter- 
esting things to observe here. 
Firstly, there is a slight sup- 
pression of the critical temper- 
ature in the region of easy-axis 
uniaxial single-ion anisotropics 
{D > 0) for any E ^ 0. This 
unusual feature can be inter- 
preted in terms of strengthen- 
ing of the quantum fiuctuations 
caused by the non-zero biaxial 
single-ion anisotropy. Actually, 
the greater the ratio E/D, the 
more impressive reduction of 
can be observed. Secondly, the 
insert of figure 2 clearly demon- 
strates that the system exhibits 
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Figure 2. The reduced critical temperature 
against the uniaxial single-ion anisotropy at 
E/D = 0.0, 0.2, 0.3, 0.4, 0.5 in ascending order 
along the direction of arrows. The case E/D = 
0.0 is plotted as a dashed line for visual effect. 
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interesting reentrance in the vicinity of the critical value Dc/\J\ given by the condi- 
tion ()16|) . For the uniaxial anisotropies slightly below this critical value, one starts 
from the disordered ground state before entering the ordered phase at lower crit- 
ical temperature and finally, the ferrimagnetic long-range order disappears due to 
thermal fluctuations at the upper critical temperature. Notice that such pure reen- 
trant transition with two consecutive critical points have not been reported in the 
analogous mixed-spin system on the honeycomb lattice [|7j. This observation would 
suggest that a presence of the reentrant transition can possibly be attributed to 
higher coordination number of the spin-1/2 sites. 

The influence of the next-nearest-neighbour interaction on the phase diagrams is 
illustrated in figure 3. As one can see, the spin system remains long-range ordered for 
arbitrary strength of the ferromagnetic next-nearest-neighbour interaction (J' > 0). 
When small ferromagnetic cou- 
pling J' is taken into con- 
sideration, however, the criti- 
cal temperature versus uniax- 
ial anisotropy dependence can 
be characterized by the 5*- 
shaped curve indicating an oc- 
currence of the reentrant tran- 
sitions with three consecutive 
critical points (see the insert in 
figure 3). It is noteworthy that 
this reentrant behaviour gradu- 
ally vanishes as the ratio J' /\J\ 
strengthens. Even more inter- 
esting situation emerges when 
assuming the antiferromagnetic 
next-nearest-neighbour interac- 
tion (J' < 0). In the latter 
case, the displayed critical lines 
imply an existence of reentrant 
transitions with either two or three consecutive critical points before approaching 
Dc/\J\ given by equation (fTHj) . Notice that a range of the critical boundaries cor- 
responding to the reentrant transition with two (three) critical temperatures over- 
whelms at relatively smaller (greater) values of | J'|/| J|, while the reentrance com- 
pletely disappears from the critical lines above certain boundary value of | J'|/| J|. 
Finally, it is worth mentioning that similar influence of the next-nearest-neighbour 
interaction on the criticality is observable even if a non-zero strength of the biaxial 
single-ion anisotropy would be considered. 

3.3. Temperature dependence of magnetization 

In this part, the temperature dependence of magnetization (to be further abbre- 
viated as m(T) curves) will be explored in particular. When both kinds of atoms 




-1.0 -0.5 0.0 

DI\J\ 

Figure 3. The reduced critical temperature 
against the uniaxial single-ion anisotropy at 
J'/\J\ = -0.2, -0.1, -0.05, -0.025, 0.0, 0.025, 0.05, 
0.1, 0.2 in ascending order along the direction of 
arrows. The case J' /\J\ = 0.0 is plotted as a 
dashed line for visual effect. A dotted line shows 
the line of compensation temperatures. 
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(a)£/£) = 0.0 (b)£/D = 0.2 




Tl [J\ k^ Tl \r\ 



Figure 4. Temperature variation of the resultant magnetization as the uniaxial 
single-ion anisotropy varies and the ratio J' /\J\ = 0.0. The parameter of biaxial 
single-ion anisotropy is fixed to be: (a) E/D = 0.0; (b) E/D = 0.2. 



are stochiometrically taken into account, the resultant magnetization reduced per 
one lattice site of the mixed-spin diced lattice can be defined as m = + 2ms |/3. 
Figure 4(a) depicts the effect of uniaxial single-ion anisotropy on the shape of m(T) 
curves. The downward curvature observed for strong negative uniaxial anisotropies 
can obviously be attributed to a more rapid thermal variation of the sub-lattice 
magnetization m^, which is thermally more easily disturbed than tua under the 
condition < 0. In accordance with the phase diagrams displayed earlier, one 
finds m{T) curves exhibiting the reentrant transitions with two consecutive criti- 
cal points slightly below the critical value Dc/\J\ given by the condition (fTT)|l . For 
comparison, figure 4(b) illustrates how the biaxial single-ion anisotropy affects a 
temperature variation of the resultant magnetization when the ratio E/D = 0.2 is 
fixed and D/\J\ varies. It is quite evident that the most ultimate difference between 
the m(T) curves depicted in figures 4(a) and 4(b), respectively, rests in the zero- 
temperature limit of the latter ones that do not start from their saturation value 
owing to the quantum reduction of the sub-lattice magnetization ms caused by the 
biaxial anisotropy. 

Next, the effect of next-nearest-neighbour interaction J' upon the shape of m{T) 
curves will be clarified. The effect of biaxial anisotropy is neglected here for sim- 
plicity, however, it worthy to note that the curves displayed below remain qualita- 
tively unchanged even if the biaxial anisotropy term would be non-zero. Figure 5(a) 
shows typical m(T) curves under the assumption of small ferromagnetic interaction 
J'/|J| = 0.05. The most remarkable differences compared to the m(T) curves dis- 
cussed earlier can be found in the vicinity of the boundary value Dc/\J\ = —1.5, 
which is strong enough to force all the spin-1 atoms towards their non-magnetic state 
S' = in the ground state. For illustration, figure 5(b) depicts the sub-lattice magne- 
tization of two most striking dependences in addition to the resultant magnetization. 
As one can see, there appears m(T) curve with one compensation point if the uniaxial 
anisotropy is chosen slightly below its boundary value (for instance D/\J\ = —1.6). 
The compensation phenomenon arises on account of thermally induced increase of 
the sub-lattice magnetization niB, which overwhelms the sub-lattice magnetization 
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(a) J'I\J\ =0.05; EID = Q.Q 



(b) J7 |J1 = 0.05; EID = Q.O 
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Figure 5. The effect of small ferromagnetic next-nearest-neighbour interaction 
J' /\J\ = 0.05 on the temperature dependence of magnetization. Figure 5(a) shows 
the variation of the resultant magnetization, Figure 5(b) depicts in addition to the 
total magnetization (solid lines) also the temperature variation of the sub-lattice 
magnetization (dotted lines) and (dashed lines), respectively. 
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Figure 6. The same as in Figure 5 but for the case of small antiferromagnetic 
next-nearest-neighbour interaction J' /\J\ = —0.05. 



niA above the compensation temperature. When the uniaxial anisotropy is taken 
to be even more negative, there emerges m(T) curve with three consecutive critical 
points instead of the one with the compensation point (the case D/\J\ = —1.63). 
Notice that further decrease of the uniaxial single-ion anisotropy shrinks the tem- 
perature range between the second and third critical point until the ordered phase 
completely vanishes from this higher-temperature region. 

At last, we shall briefly discuss the effect of small antiferromagnetic next-nearest- 
neighbour interaction {J'/\J\ = —0.05) on the temperature dependence of the mag- 
netization as shown in figure 6. According to equation (fTHj) . the critical line ter- 
minates due to the spin frustration already at Dc/\J\ = —1.475. It is quite ap- 
parent from figures 6(a) and 6(b) that the magnetization exhibits two reentrant 
transitions for uniaxial anisotropics slightly below this critical value, for example at 
D/\J\ = —1.5. Contrary to this, the resultant magnetization may show the reentrant 
transition with three critical temperatures when the uniaxial anisotropy is selected 
slightly above this boundary value (the case D/\J\ = —1.47). Whether the reentrant 



10 



Reentrant transitions of a mixed-spin Ising model 



transition with two or three critical points is observed, however, depends very sensi- 
tively on the relative magnitude of all interaction parameters. It appears worthwhile 
to remark that both kinds of reentrant transitions disappear as the magnitude of 
next-nearest-neighbour interaction is considered to be strong enough. 

4. Concluding Remarks 

In the present article, the mixed spin- 1/2 and spin-1 Ising model on the diced 
lattice has exactly been treated within the generalized star-triangle mapping trans- 
formation. The main focus of the present work has been aimed at the influence of 
the next-nearest-neighbour interaction, uniaxial and biaxial single-ion anisotropics 
on the magnetic behaviour of the system under investigation. It has been found that 
the single-ion anisotropy may possibly cause the reentrant transition with two criti- 
cal temperatures only. On the other hand, the ferromagnetic next-nearest-neighbour 
interaction may lead just to the reentrant transition with three consecutive critical 
points, whereas the antiferromagnetic next-nearest-neighbour coupling may be con- 
sidered as possible cause of both kinds of reentrant transitions. 

Another interesting finding to emerge here consists in the fact that the biax- 
ial singlc-ion anisotropy has been confirmed to be a possible source of the quan- 
tum reduction of magnetization. Our further observations suggest that the biaxial 
anisotropy cannot qualitatively change the character of reentrant transition caused 
by the uniaxial single-ion anisotropy term or the next-nearest-neighbour coupling. 
Finally, we should remark that a combined effect of uniaxial and biaxial single-ion 
anisotropics with the next-nearest-neighbour interaction is responsible for manifold 
temperature dependences of the resultant magnetization by assuming the ferrimag- 
netically ordered spin system. 
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